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ABSTRACT 


Predictive model was developed to detect insulation failure and pipeline leakage in natural gas transmission line using simulation 
software known as EXTENDSIM9 SUIT. The computational fluid dynamics parameters for the flow of compressible gas in a pipeline 
includes density, temperature viscosity, pressure etc which changes along the transmission line due to frictional resistance, and may 
lead to fogging/ hydrate formation or sudden degradation of pipe materials like erosion which gradually leads to corrosion or rust, 
pitting etc. This anomaly is one of the major contributions to pipeline failure thus this developed model can be used to monitor and 
predict insulation failure along the natural gas transmission pipeline since there is a lot of the discrepancy between the 
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Computational Fluid Dynamics (CFD) parameters and their arrival components. Sensors where used to read and monitor fluctuation 


of fluid flow dynamics properties which was set to alert the indicator otherwise known as the decision block to analyze the results or 
records obtained from the constraint variable blocks which reads 12.5Km due east as the failed distance. From the sensitivity 
analysis result we observed that the model can predict insulation failure along the transmission line with 92.2% accuracy as it is 
sensitive to both small and large parameter difference and from our findings, the temperature increased exponentially along the 
failed area. When these anomalies are detected by this model; it triggers a maintenance shutdown on the affected area and record 
“time between failure and time between repair” by using (FDIMs) Fault Detection/Isolation and Maintenance Shutdown). Thus this 
model can be utilized to approve the readings of an Advanced Fiber-optic Distributed Temperature Sensors (DTS) or Distributed 
Acoustic detecting (DAS) in a gas transmission line. 


Key words: Predictive model, detect, insulation failure, pipe leakage, natural gas transmission pipeline, simulation software 


1. INTRODUCTION 


Natural gas is a naturally occurring hydrocarbon gas mixture which is combustible as well as the main sources of energy for 
domestic and industrial uses. The composition made of mixtures of hydrocarbon, such as ethane, propane, butane, and pentanes as 
well as the presence of water vapor, hydrogen sulfide (H2S), carbon dioxide, helium, nitrogen, and other compounds. The natural gas 
is classified as fossil fuel which is utilized as a vital source of energy for heating, cooking, and electricity generation. It is additionally 
utilized as a fuel for vehicles just as a chemical compound feedstock in the manufacturing of plastics and other monetarily 
imperative natural synthetic substances. Research conducted revealed that fossil fuel based natural gas is a non-sustainable 
power/energy source, because its end product after compete combustion are carbondioxide, water, heat etc [1]. 

Investigation conducted on transmission pipes revealed the significance of pipe diameter, the few characteristics and behavior as 
well as the effect on the functional parameters of transmission principle guiding the control process [1]. Further considered the 
significance of pipe dimension effect as well as material selection in pipeline construction in control system for detection of 
insulated failure which leads to leakages in pipeline; other factors that contribute to the pipeline failure include hydrate formation 
formed as a result of change in temperature in the internal flow of fluid in the pipeline [2]. The change will influence the pressure of 
flowing liquefied natural gas in the pipeline as well as the viscosity, density and other functional parameters will change as well. 

Research work conducted on the significance of coating a pipeline reveals the prevention of corrosion and resting as a result of 
contract of external surface of the metal with the presence of moisture content, microbial activities and other environmental factors. 
Initially, major coating materials used on pipes are coal tar enamel but in the recent time the major materials used are fusion bon 
epoxy and other polymer materials [3]. 

The word insulation” can be described as a process that enhances protection or isolation of an existing system. The concepts of 
protection or isolation in pipeline are found necessary because of continuous failures currently observed by the oil companies. The 
insulation process of their pipelines enhance safety quality in terms of detection of regions, which is likely to fail as a result of 
malfunction in the process fluid, which in turn influence the integrity of the pipeline material [5-6] 

The new technology adopted by the various coating company using polymer material product has rendered the initial approach 
of using tape coating system and coal-tax enamel system, because of the deficiency found in terms of insulation as well as the 
constrain when large diameter pipe is to be insulated. Research conducted on using tape coating and coal-tax enamel is more 
effective in small diameter pipe as well as fail when the environment is wet compared with coating using polymer materials [7]. The 
polymer materials possess high resistance to wet environment as well as reduce the failure rate of the pipe and enhance 
performance of the system. Research conducted by various Scientifics group revealed that the application of coal-tar poses high 
level of health challenges to the environment [8]. 

In natural gas transmission pipelines, insulation is performed in other to reduce heat loss from the pipe. However, some of the 
major reasons of insulation failure includes: Natural decay or deterioration because of maturing, which is quickened by excessive 
heat, moisture and dirt, chemical deterioration, mechanical damage and sunlight are main causes of insulation failure. Other causes 
may include excavation, incorrect installation and operation, pipeline vandalization, accidents etc. 

This is widely used in chemical engineering to detect failure in chemical plant, petrochemical as well as in information 
technology, disaster recovery, meteorology, capacity and city planning, etc. (Schiff, 2002) however, it is found useful in detecting 
insulation failure in terms of pipeline leakage, which enhance reduction pressure, temperature as well as hydrate formation in 
pipeline. It consists of various factors which influence the future behaviors or results of the future failed system using extendsim9 
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software. This is achieved mostly by or relies on mathematical model simulation model built on the software as well as the system to 


identify the discrepancy between (functioning system) and faulty part of the system [8-9]. 

Failure of an insulator in natural gas transmission pipeline causes pressure drop along the pipeline. This is due to the fact that 
heat has been transferred either by conduction or convection from the liquid to the wall or from the encompassing to the liquid in 
this way causing operational instability which affects the flow and gradually leads to hydrate formation, which causes blockage along 
the pipeline. However, the following characteristics will be observed such as change in pressure, temperature and density which 
causes erosion, gradually leads to the exposure of corrosive (H2S) chemical substances to the pipe material and when these 
anomalies are not detected on time it exposes the pipe metal to corrosion, pitting, pipe failure and eventually explosion. The aim of 
this research work is to use existing simulation software named Extendsim9 to detect insulation failure and pipe leakage in natural 


gas transmission pipeline. 


2. MATERIALS AND METHODS 

Extendsim9 Simulation Model 

Just like other simulation technology, Extendsim9 is a multi-purpose simulation software with combination of mathematics, 
computer, physics and engineering, and graphical user interface using a combined logical representation to select few resources put 
together to build entire systems to a particular function [10]. To solve real — time pipeline flow on fluid mechanics and dynamics 
problems, it is ideal to have good knowledge of the fluid flow behavior as well as its properties to arrive at a better solution. 


Table 1 LNG pipeline layers (NNPC) 


Pipeline layers 


Layer Material Wall Heat Thermal Density 
Thickness capacity conductivity kg/m? 
mm j/kg.k w/m.k 

Transmission pipeline Steel 12,7 500 50 7850 

Corrosion allowance Steel 1.5 500 50 7850 

Anti-corrosion coating Polyethylene 3.2 1020 0.02 1293 

Insulation Polyurethane 50 2300 0.03 950 

Air 10 1500 0.025 32 

Outer shell Steel 12.7 500 50 7850 

Transmission Pipeline external Diameter 648.6mm 

Concrete thickness 50 

Concrete thermal conductivity W/mK 2.79 

Ambient temperature °C 15 


Mathematical Model Development 

In developing a model, there are certain parameters to identify the type of flow before development; 

° Compressibility: Compressible flow otherwise known as variable density flow, this is as result of dual nature of gases as the 
density, temperature changes, depending on the operational condition. For air, 10-5 m2/N at 1atm. This is given by 
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Ko = a ae in m*/N 


p\dp (1) 


° Prandtl number: This is used to determine the relative viscosity of the flow which satisfy it viscous or inviscid flow. For 
typical gas flow it values ranges from 0.7-1 


2165x1.9534x107° _ 


1.4 (2) 
0.030 


Pr = Cpp/K => 


Where, K = thermal conductivity, 1. = dynamic viscosity, Cp = specific heat capacity at constant pressure. 
Knudsen number is the ratio of fluid mean free path to the macroscopic length scale of a physical system. 


Kn = ZL (3) 


Where, X = ———==— (4) 


fs) 
oe ap (5) 


Ox 


Where, X= mean free path, p = dynamic viscosity, R = universal gas constant, T = temperature, L = length (characteristic length) 
Hence, If Ky < 0.01 continuum will not hold. 


Beyond this crucial range or Knudsen, number; If 0.01< Kn < 0.1 = slip flow, 0.1 < Kn < 10 = transition flow, Kn > 10 = free molecule 
flow. From our calculations, K, = 0.7. Therefore it is a continuum flow. From these; we can deduce that the flow is continuum, 
compressible and a viscous flow. 


Modeling Considerations 

This research considers continuum, compressible viscous flow along the natural gas transmission pipeline occurring in NNPC. This 
work tends to develop flow models that can be used to detect insulation failure and pipe leakage on the gas transmission line. The 
predictive models will be derived from the conservative continuity, momentum and energy equations, fluid/aerodynamics and heat 
transfer equations. 


Mathematical Model Assumptions 

The following model assumptions were considered during the research work as stated below: unsteady state because of the faulty 
part of the system; we intend to achieve steady state, gas flows are transonic or turbulent, heat loss by convection from the inside 
pipe flow, conservation of mass, momentum and energy, heat loss by convection to the outside pipe, heat conduction across the 
insulator, heat conduction across pipeline thickness, thermal resistance between the fluid and the insulated wall, the wall itself and 
the resistance between the wall and the insulator. 


Basic Pipeline Equations Used For Modeling Gas Transmission Pipeline Flow 

Fluid mechanics and aerodynamics offers a systematic structure that underlies physical disciplines which embrace empirical/ semi- 
empirical laws derived from flow measurement which is used to solve practical problems. The solution to a fluid dynamics problem 
typically involves the calculation of various properties of fluid such as flow velocity, pressure, density and temperature as a function 
of space and time which describe the behavior of fluid being transported at all points within the pipeline flow system. The basic 
equations are the continuity equations, momentum equations and energy equation [11]. 
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Continuity Equation Development 


The Continuum Equation: This implements the law of conservation of mass in a pipeline system which means that the actual mass 
flow in and out of the pipeline section is equal to the rate of change of mass within that section. 


pl 


x 


p2 
@) to 


Z 


Figure 1 3- D Cartesian coordinate flow 


Consider a flow in 3-dimensional Cartesian coordinate, the mass balance becomes 


Rate of change of \ {Rate of inflow Rate of outflow 
( mass within the of mass into }=| of mass from (6) 
pipeline the pipeline the pipeline 


Analysis of equation (6) mathematically and considering total density flow from point (1) as initial to point (2) as the final flow. 


Pi Ce eae, (7a) 
P2 (x. Yo» Zreto) (7b) 


Using Taylor's series to obtain the initial and first derivatives of the flow, and also truncating higher derivatives gives: 


fle y,z,t) _ faen VorZasta) + ae AG, Yor Zrto) 90) Cee iets) 


Fis ys zeihs) + (t-1,) (8) 


Applying equation (8) into the equation (6) yields: 


0p 0p op 
Pr = pt—(x, —4,)+—(.- 9, + - 4) 
ox, oy, 0, 
(9a) 
1p 
+—(t =i) 
ty 
Dividing equation (9a) by tz — ti throughout; we have 
Pe Pi 2 [se }=a), op banm),,(@){a=a) , (9b) 
ee 2 Ox) in =1 Cy} 15=h Ce) 13—1, Ot 
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Taking the limits tz — t; on the RHS of equation (9b), we have 


P2— Pi = i 
a 7 At 
Lim t-t 


t, —ty 


X, —X, = Lo Tk 
Where; 2 u = V. Y27 = V, 2 -_ ¥, are the velocities of the different Cartesian 
t,—1, t,—t, “(by -h 


coordinates in x, y and z-directions respectively. 


Substituting these ( oh : Vi; V,, V.) into equation (9b) we have, 


oe = V OP 59, Pay oe OB 


At Ox Oy “oO Oat (10) 


where wie = the time rate of change of density of a fluid as it moves through a space. Equation (3.10) can be written as 
A 0 

== V yy, +1V, Y% 4v oA (11) 

At Ot * SOx y Oy 2 /OZ 


Considering the different directions of flow and substituting the flow velocity vector into equation (11) we have; 


A Z Oe, o ek Fee 
re = apt Ve Yay IV, Yagt kV, Y 
oe = wa +07) 
Where, V = PV, +jV,+kV,) =V = EV, +jV,+KV, 


Is the flow velocity vector and 


> 


Val f+ I Vy on k Vp, Is the divergence vector 


a 
Where, us = local derivative, which is physically the time rate of change at a fixed point and [vv = convective derivative or 


divergence velocity which is physically the rate of change due to the movement of the fluid element from one location to another in 
the flow field where the flow properties are spatially different. Or it is the time rate of change of the volume of a moving fluid 
element per unit volume. Within the control volume, the total change of density in a close loop is equal to zero. This means that: 


Ap _ 
At 
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Hence, equation (12) becomes: 


a2 + vf fe) v| = 0 
Ot 


(13) 


Thus, equation (13) is conservation of mass for continuity equation flow. 


Momentum Equation Development 

The momentum equation: This implements Newton's 2™ law when applied on a moving fluid element; which says that the net 
force on the fluid element equals its mass times the acceleration of the fluid element. Since it’s a vector relation which can be split 
into x, y, and z-direction. The two major type of forces described by Newton's second law: Volume forces - This has to do with 
gravitational force which acts directly on the volumetric mass of the fluid element. Surface forces - This has to do with pressure 
distribution acting on the surface which is imposed by the outside fluid surrounding of the fluid element, and inner frictional forces 
or viscosity force acting on the surface through pushing or tugging. 


(x+ Ax, y + Ay, z+ Az) (PV) c+a2 


(2V,). 


(x, y, Z) 


Figure 2 Momentum Balance for Elemental Volume on a 3-dimensional Cartesian Coordinates 


For Momentum Balance (Momentum Transfer) of Fluid Flow in Pipeline; 


Rate of momentum Raeotimonentan Rate of outflow of 


Accumulated within : he ~ momentum from 
Entering pipeline 


Pipeline Pipeline 


Summation of forces 
+ acting on the 
pipeline system 


(14) 


O 
Ax A — 
y Aza (pv,) 


Rate of accumulation = (15) 
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Ax Az((pv2), )+ Ax Ac(p(v, v,),) 


Rate of inflow of momentum = (16) 


_ Ax Ay((pv,),) 


Rate of outflow of momentum = 


Ax Az((pv,.V,) car Ax Az(pv, v,) + Ax Az(pvz VZ) rast Ax Ay((pv.).) (17) 


yrAy 
Sum of forces acting on the system (the stress, pressure forces and the gravitational forces acting on the fluid) 


= (c.., ), Ay Ax + (p, =a )Ay Az + pgx sin 0 Ax Ay Az (18) 


Le =(60 +p. +F. — Summation of the forces acting on the system 
L x P. y g y 


Combining these equations into equation (14) gives; 


a 
Ay Axe? = —Ay Az O(pv ,.v,)- Ax Ay O(pv,.v,) 


Ax Az O(pv,.v,)- Az Ay (z,,) + pg Sin 9 - Ay Az (Pyne — Px) = 19) 


Ax Ay Az 


Dividing all through by gives; 


0 
do) (An) Am), Ao) ia 2 2a) Bg Sa 
at a. @ a. | a, 


x y 


—2 
Recall that P8,=Pfi X= 
a (21) 


Applying Navier-Stokes rule on equation (20) Conservation of momentum equation and incorporating all the forces acting on the 


system yields: 


2 > 
alr.) alow ), Alor. low.) = ® : [an +2uU aI ce | 
Ot ox oy oz Ox Ox ox } Oy| \ ax dy 
0| Ov, | Ov, (22a) 
aE 7 +, 
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a( »)  Alev.r,) Alor?) | aov,v.) = Op ats ce} Le os } 
ot Ox oy Oz Oy Ox Ox = oy 
> OV fs) ov, Ov, (22b) 
Cf | AVV +2— : » 
+%| ae Oy aul Oy . Oz JJ 
a(ev.)  a(ev.v.)  AeV,V.)  alov2) _ ap , 0 ( av. | OV, ) 
ot Ox Oy Oz Oz Ox Oz Ox 
ov. OV, (o) > OV. (22c) 
© : a AV.V +2 < | of. 
“Yo af Oy GE 2 ce ar 


Where, A = bulk viscosity coefficient, py = dynamic viscosity 


Energy Equation Development 
The energy equation: This implements the principle of first law of thermodynamics which states that total energy of a system and 


its surroundings remains constants when applied on the moving fluid element, the rate of change of energy inside the fluid element 


is equal to the net flux of heat into the element with additional rate of work done on the element due to body and surface forces. 


The Rate 


of change Net flux Rate of work done on 
of energy inside} = of heat the fluid element due 
the into the to body and surface 
fluid element fluid element forces (23) 


E + BWA) xin 
Ox 


+ Avr.) Ta) ax 
Ox 


, E + oi ax bd 
Oz 


ze vr. + Aree) Jandy 
oz 


Figure 3 Energy Balance for Elemental Volume of Fluid Flow on a 3-D Cartesian Coordinates (x-direction shown only to avoid 
clustering) 


From Figure 3, considering the work done only in the x-direction, we have: 


The rate of work done 
by the body force acting 
onthe fluid element moving 
at a velocity (24a) 


> (B) = pf.tdxdydz 


Considering work done by the body forces are due to pressure in the x-direction, yields 
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— dxdydz 


(24b) 
Also, the rate of work done by shear stress in the x-direction, we have: 


A(vatyx) dxdydz 
dy y 
(24c) 


Then the rate of work done on the moving fluid element due to both surface forces and volume forces in the x-direction is: 


[- O(v,P) 4 O(VxTxx) 4 A(vxtyx) 


O(VxT zx) 
dx dx dy + dz dxdy dz 


(24d) 


Therefore the total rate of work done on the moving fluid element is the combination of surface forces and body forces in the x, y, 


and z-direction. 
From equation (23), we have 


The net rate of work done 


on the fluid element = [- a ae A(vyP) 4 es he O(VxTxx) 4 O(VxTyx) 4: O(VxTzx) gfe A(vyTxy) 4 A(vytyy) 4 A(vytzy) 4. 
dx dy dz dx dy dz dx dy dz 
due to body and surface forces 
O(vzT xz) A(v,Tyz) O(vzT zz) > 
arte nae +2Cstee)| dxdydz + pf.Vdxdydz (24e) 
. A@yP) , A(vyP) | AwsP) _ = 
Note that: — oe a ae V.(Pv) (24f) 


Also from equation (23), the two reasons for heat flux into the fluid element is as a result of; 


i.Volumetric heat e.g. absorption or emission of heat 


ii. Heat transfer across the surface due to temperature gradients, (conduction) 


The mass of the 
Nou uat { moving fluid ai) = ponayiod 
For volumetric heat . 
oe! o the fluid aes =r Bae (249) 


{ Heat transferred across } _ ( 


Ode 4 Ady 5 O4e 
the surface (conduction) , dy - a dxayaz (24h) 


dx 


Recall that heat transfer by conduction is proportional to the local temperature gradient 


Combining equation (24g and 24h), we have 


, (Ade 4 Ay 4 Ode 
[od (2 + ey +) dxdydz| (25) 
OT OT 4 OT 
Where,q, = ke; y= Na de ka 
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q = the rate of volumetric heat addition per unit mass, k = thermal conductivity 


Therefore, 


The net flux of a i da 4 


netateie | PUt EOD + S65) +R(ez)esiee co 


Also from equation (23), we have 


of energy inside amoving fluid energy per energy per 


the time rate of change the total energy of its internal its internal 
the fluid element per unit mass unit mass unit mass 


= Total energy = (e + =) (27a) 


When this is applied to a fluid passing through a fixed control volume and the mass of the moving fluid = pdxdydz, we have; 


the total energy of . 
the moving fluid | = p—(e + -) dxdydz (27b) 
per unit mass 


Combining equations (27b), (26) and (24e), we have 
Energy equation becomes 


pa(et+>) =[pg2(k2) +2 (c2) +2 (c2)] -?- BlvyP) _ BveP) | Aorta) 5 A%ytyn) 5 Bester) 5 Weytey) , Aeytyy) 


dx oy dy az dz Ox oy Oz Ox oy Oz Ox oy 
A(vytzx) O(VzT xz) O(VzT xz) A(ztyz) O(2Tzz) ro 
D 28 
Oz Ox Ox si oy + Oz - pP y ( ) 


Note that equation (28) is in non-conservation form 


Transferring equation (28) to conservation form in terms of total energy, we have: 


From the LHS of equation (28), 


a(e+>) a(e+=) v? 
p+ = p++ pi.v(e+7) (29a) 
But, 
v2 v? 
plore) _ elo) (4 ye (29b) 


From flow vector velocity regarding divergence 


Conservation of energy 


2 2s 
of ex Poul ee aie = piv S(rZ). 2 or +2(eZ) OV. p 
ot 2 2 Ox ax) dyl dy) Ol &) a& 
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AV, p)_alv.p) , air aaa its 7 Va mae maae F aV.t,, 


XXX 


oy Oz Ox dy Oz Ox dy Oz Ox 

ov. 30 

gee Mite, piv (30a) 
dy OZ 


Equation (30a) is the conservation of energy and can be re-written as 
= 


+ VEG = S-V.g + pf.b+V.ty6 (30b) 


Where e = internal energy, t = shear stress, Where Q = total flux of heat into the element 


> 5 
J .»f = body force per unit mass acting on the fluid element in x-z direction 


Application of Continuity, Momentum and Energy Equations on Pipeline For 3-dimensional cylindrical coordinate flow 


Where, 

X = rcos@ >r 
Y= rsin@ —6 
Z= zZ —Z 
8 = _ Arctan (y/x) 


Velocities for 3-dimensional cylindrical coordinates, 


Vx = V; 
Vy = Vo 
Vz = Vz 


For enthalpies in 3D cylindrical coordinates considering the r-direction as reference coordinates; 


(31a) 


(31b) 
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hy = il he = r, h3 = | (31c) 


Converting the Conservation of continuity equation (13) into 3-dimensional cylindrical coordinates; 


Continuity equation becomes => 


ILIAC @ : 
at - a eve) A, (eV.) 0 (32) 


Also, Transforming Conservation of momentum equation (24g), (24h) and (241) into 3-dimensional cylindrical coordinates; 


Momentum equation becomes =r - component 


+ pgr 


V V V V y ; ov Op 10 1 Or. T OT 
2 ray 7) # 6 6) - 6 V ’ ) r@ 06 p 
(33a) 


oO 0 'or ro0oé r ‘& orror ”" r0@ r & 


8 -component 


10 loys 
pear ee Mey Me) 7 ae 2 (25) He pp 


a? a 6]: a nn dOr 00 r or dz (33b) 
z-component 
j ov. y ov, Yo ov, 4. ov. \__ p10 Gg ig es 
Ot or r 00 * & Ozror' ~ rod & (330) 


Transforming Conservation of Energy equation (30a) into 3-dimensional cylindrical coordinates, we have; 


OE, 60 
ot Ot 
(EV, + pV, -V,T,6 —ViT 99 —V 06, +4y)+%, (EV, + PV, =V iT Vibe, Vie +q, )=04) 


0 
— ALLY, + fo + FY. )+— (EV, + BV, Vit Voto Vita +4, }+ ee 


Where; 
e+V? 
E, = 7 pe 
(35a) 
But, 
Ppf=ps8 
(35b) 


Where stress tensor 1 
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(35c) 


Energy dissipation function ® 
2 
aVv,\' 4/a, V av, \’ a (av av.) (av, av.) 
cy) = 2 a (en ec ee re ae iY — | 4 j—— 
4(%) ee 7 (=) AZ r ye) fe ) 


év, 1éV.) 2 \ i) 
. TV 
if ar ) cen 


(35d) 


Enthalpy => eee ane: 


- (35e) 


Combining equation (35a) — (35e) and substituting in equation (34) we have; 


Ot r or r 00 Oz 


eV, 16V, V, eV. a(Vv,) 1éV, eV. OV, 
T, ——+T gy | -——+ +7, —|]-|7,,] 7 + +7,.| + 
or ro@g er Oz or\ r r 00 or Oz 


Sle k))+2 5 (oany,)+4< (ohive)+ 2 (ohv.)+{ 2S : 
Yr or 
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10V, OV, 1 1 Or T OT 1 Or 
rs 2 ge ey 20 te ee ee oy | 
re(* 00 ] (F%br r 00 r @& [3 00 ) 


1 OT 1 1) 1é@t,.. at 
a 5) 2 & -0 & We 
+ 2 / or (- r,)+ a v{2 ole, += |+2 s + a = Energysource 


(36) 


Based on the above differential equations the internal flow parameters are obtained using software. 


ra) + k > > > > 
Recall that, aple +k) + v{ ol + k)V+ pV-t,.V+ ‘| = 0 
ot (37a) 
Using continuum equation, if only kinetic and internal energy are considered 
E 
{"%)} 
OE > 
p22 a Py any 
At Ot (37b) 
E 
AY . 
; %) : jae, av) 
At At At ais 
But from velocity vector V we deduce 
A 7 > 2° OS > > 
ee = pf V-Vp.V+(V.z,)V 
At (37d) 
Combining equation (37b), (37c) & (37d) and substituting in equation (30b) above we have; 
A > ra] > > 
ps fv) = “Maggs vr, ¥ -(V.z,)V 
At Ot (37e) 
But h = vx, ¥ — (V-z,) 
and (37f) 
> ra] ; 
V = T, = 
oy (379) 


Substituting equation (3.37f) into (3.37e) we have 


p+ [vv] = —-V.q+@® 
t (38a) 
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Finally, equation (3.38a) can be written 


Ah Ap oO 
ge =, Pe gee 
At At Ot (38b) 
Equation (38b) is the conservation of energy 
> > 
f = force per unit mass, g = body force (acceleration due to gravity vector, e = internal energy per unit mass, E, = total energy 


OEt 
per unit mass in the control volume, —— = the rate of increase of total energy in the control volume, V.EV = the rate of total 


Ot 
y 
energy lost by convection per unit volume, —— = the rate of heat produced per unit volume by external force, o f .V = work 
t 
done on the control volume per unit volume by the body force, V.g = the rate of heat lost by conduction per unit volume 
> 
through the control surface, Vij Vs = work done on the control volume per unit volume by the surface force 
Heat transfer in gas transmission pipelines 
A A 


Figure 4 A simple insulated pipe (well- insulated pipe) 


The process involves the described functional parameters as stated below: the heat loss by convection from the inside pipe fluid, 
heat loss by conduction across the pipe thickness rz —1r1, heat loss by conduction across the insulation thickness r3-r2, and heat loss 
by convection to the outside pipe fluid 


Figure 5 Systematic size view of the cylindrical insulated pipeline 
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From Fourier’s law, heat transfer by conduction is given by; 


_ — kot 
oe Ox 


(39a) 


For conservation of energy on 3-dimensional cylindrical coordinate equation (39a) becomes 


— kot k 
q..= 4 eee q. oes! 
Ox Oy Oz 


(39b) 
Where, k = thermal conductivity 
Recall that; 
\ 
a aq, dq, 94, @ 
q.-|9,+ tar loa = =- “a... jg=| 2a eo. 
Ox Ox , ox Oy) Cz " ” 390) 
Recall also that 
_ ar 
1 = Sy 
(39d) 
Where, AT = T11- Te (inlet temp-arrival temp), 7; = thermal resistances 


Thermal resistance (R) is defined as the ratio of change in temperature to the associated rate of heat transfer: the heat transfer 
between the fluid and the wall, the insulated wall itself is one resistance and the heat transfer between the insulated wall and the 


outside fluid 


But the overall heat transfer coefficient can be calculated as 


RA 
(40a) 
recall thatR = Vy 
(40b) 
1 1 (40c) 
= c 
UA hei Ai t+ U(Sp/knAn) + 1/AeoAo 


Where, R = heat transfer resistance m*K/W, U = the Overall heat transfer coefficient W/m? K, K = the thermal conductivity of 
material in layer (W/(m.K)) Sn, = thickness of the pipe mm, h¢io= inside or outside wall individual fluid convection heat transfer 


coefficient W/m? K 
But total heat loss from the pipe 


Q = MCpdT > Q = MCp(T;-Te) 
(41a) 


Where, 
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= A, 8 
Tg = T; ee 
(41b) 
AT =T, —-T,, 
(41c) 
T. -T 
T, = i E 
2 
(41d) 
Where, 1; = inlettemperature, Te = exit/arrival fluid temperature 


Combining equation (40c) & (41c) and substituting in equation (39d) above; we have, 


\ ee § 


i E 


r, r, r, 
l eo ee l 
Qahr.L° 2nk,L° 2nk,L nk, L 2ah,r,L 


q= 


Equation (42a) is written in C++ programming language Code used to simulate the total heat loss of the normal insulated system 


(without failure) 


q= T; ~~ T, 
Be ” 
l +L A +L a + : 


Qa  °2akL * 2ak.L QahrL 


Thus Equation (42b) is written in C++ programming language Code used to simulate the total heat loss of the failed insulated 


system. 


Predictive Model Development 
For Pressure 
Consider the faulty part of the flow system to be one dimensional in the r-direction, neglecting other flow directions 


This model is from Equation (33a), we have 


ov. Wve Vi OY ap 10 LOrs. ty OF . 
fe) = — po 4VY—* =— rT, )+ —-—S+—=+ pgr Therefore, equation 
Ot or r0oéO fr * OZ Or r Or r 00 r_ © 


(33a) is reduced to; 


OV OV, 1 
p—t+V,— =- De (rc, )+ psr (43) 
Or Or r Or 


Assuming at steady state, the gravitational (body and volume forces) and surface forces are neglected, equation (43) becomes; 
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OV, 
AY, ; =. 2 (44a) 


udu dP udu dP 


wr +S=0 = p= -S (44b) 
Hence equations (44a) and (44b) are similar 
From equation (44a) p (v, or) = -2 = pi + - =0 (44c) 
2 (P + pv?) =0 (44d) 
Upon integration; 
fH + pv?) =0 (44e) 
P+ pve =k (44f) 
At r= 0, P= PB, % =0,R+0=k,k =P, 
Upon applying the boundary conditions, Equation (44f) becomes 
P+ pu? =P (44g) 
P—P,=-—pv? (44h) 
Therefore equation (41h) becomes; 
P=P,—pv2 (45) 


Where, P, = initial pressure, pv? = compressible dynamic pressure or impact pressure, P = final pressure 


For Temperature 


This model is from Equation (30b), we have —_ + V.E,0 = - —V.qtpf.o+V. TjjV 


Assuming at steady state, the gravitational (body and volume forces) and surface forces are neglected, equation (30b) becomes; 


0 =-V.q (46a) 
Recall that from equation (11e), one dimensional flow equation (r-direction) becomes; 
30 
Veie (46b) 
Substituting equation (46b) into equation (46a) becomes; 


0 = —-V.q, 


i) 
V.gr = I=. 4r (47a) 
Where, V = divergence, q, = heat flux (rate of heat transferred through the surface per unit time J/s), 7 = flow vector 


But Qr = KA— 
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Equation (44a) becomes; 


a a, ar 
wit ska. = 9 oe 
: a ar 
Upon integration fala +ko|=0 (47d) 
ar(r —%)) +k(T —T)) = (47e) 


Where, C = constant 
At these boundary conditions: r=0,T = To, qr =Qr, To tO =C, C=Ty 


Substituting these boundary conditions in equation (47e) we have; 


Gr (7 — 1) = k(T — To) (47f) 
© (r — 19) = k(T — To) (47g) 
Therefore, T= Tp) — f(r — 1) (48) 


For Sensitivity (Temperature and viscosity) 
This model is from Equation (38b), we have 


Ah Ap Oo 
At At Ot 
Where, p = density, h = enthalpy, P = pressure, Q = volumetric heat capacity in J/m?K (heat per unit volume), < = heat 


produced per unit volume by external force J/m?K, q, = heat flux (rate of heat transferred through the surface per unit time J/s), ® 


= Dissipation function in seconds. 


Assuming the pressure change is neglected, equation (3.38b) becomes; 


An _  9Q 
Pa: = > V.q+@ (49a) 
But Vg = 44444 (49b) 
Therefore, 
Ah _ 9Q_ (9a, 44, oq 
PRE 3 (3 + 56t 2) 7 (49¢) 


Assuming there is no heat produced by external force in the system, it means that; 


0Q _ 
+ = 0 (49d) 


Therefore, p— = - (2+ “4 ot) +@ (49e) 


Considering one dimensional coordinate in r-direction, 
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Equation (49e) becomes; 


Ooh oq 


at ~ oe 


aT 
But, q= kA 


Substituting equation (50b) into the R.H.S of equation (50a) we have; 


ES} 
| 
| 
y 
S 
oD 
2| 
y 
$ 


At these boundary conditions; t = 0, T) =T, m =0, r=L 


But  Oh= CydT 


Substituting equation (50d) into equation (50c) we have; 


CpdT _—siKAOT. 


at ~ arar 


Note that for one dimensional coordinates; 


0 = u{(%2) - Eno 


Input Parameters 


Table 2 LNG operational data (NNPC) 


Operational data 


Transmission Flow rate Inlet pressure 
pipelines m?/day KPa 

1 2370 1098 

2 3120 1108 

3 3533 1113 

4 4125 1181 

5 4330 1280 


Gas properties 
Natural Gas gravity kg/m? 
Gas specific heat capacity J/kgK 
Boiling point 
Freezing point 


Flammability limits 


Inlet 


(50a) 


(50b) 


(50c) 


(50d) 


(51) 


Outlet 


temperature °C __ pressure KPa 


50.3 
51.1 
50.5 
50.4 


51.6 


652 
656 
615 
625 


402 


0.489 

2165 

-161.5°C at latm 
-182.6°C 


4-5 by volume in air 


Outlet 
temperature °C 


42.80 
44.45 
45.12 
45.85 


41.51 
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Gas density 


Ignition temperature 
Carbon content by weight 
Hydrogen content 
Oxygen content 
Hydrogen/carbon atomic ratio 
Relative density 

Octane no 

Methane no 

Viscosity u 

Universal gas constant 
Gravitational constant g 


Thermal conductivity k 


1.0926517kg/m3 
538°C 1atm 

73 

24 

0.4 

3.0-4.0 

0.72-0.8 at 15°C 
120-130 

69-99 
1.953x10°kg/m/s 
8314.4)/kmol/K 


9.81m/S? 


0.030W/m2.k 


Plate 1 Simulation Model to Detect Insulation Failure in Natural Gas Transmission Pipeline 


Plate 1 shows the prediction of insulation failure along gas transmission pipeline using a simulation model called Extendsim9 


software. 


Ter 


re | 5 6 


ARTICLE 


Chart 1 Computational flow chart 


Write CFD 
programs 


Input parameters 


Compute simulation variables e.g 
temp, pressure, etc 


Display 
results T, 
P, V etc 
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3. RESULTS AND DISCUSSION 


Overview 


Predictive model to detect insulation failure and pipe leakage in natural gas transmission pipeline has been modeled and predicted. 
Normally, whenever there is a flow of compressible gas in a pipeline; density, temperature and pressure changes along the line due 
to gravitational/surface forces which may lead to sudden degradation of pipe materials like erosion which exposes the pipe to 
corrosion or rust, pitting etc. This is one of the reasons pipes undergo insulation process before installed in natural gas transmission 
lines in other to minimize hydrate formation, accumulation of corrosive chemical or substance like H2S on the pipe material such 
that heat will neither be lost nor gained. 

However, when there is sudden failure or failure due to aging of insulation materials on these pipelines, it poses some threats on 
the operational impacts whereby heat maybe lost or gained depending where the insulation failed either internal, external insulation 
or both which may cause the gas to flow at a lower temperature which may lead to condensation and freezing of gas along the 
pipeline, thereby causing hydrates formations, which will in turn cause pipe and equipment obstructions. When these anomalies are 
not detected on time, the possible adverse effects are explosion. In addition, insulation failure with adequate effect of high 
temperature also exposes the pipe to both internal and external corrosion, pitting and leakage. 

The inability to detect insulation failures along the natural gas transmission pipelines on time has been a great challenge to the 
oil and gas industries. Although there has been series of transmission pipeline technologies which is used to monitor and detect 
leaks, as well as pipe blockage and can only detect pipe failures with minimal attention to insulation failure when pipes’ materials are 
still intact. These inspections are carried out periodically as a survey (once in every 5 years), and thus; insulation may fail which may 
cause operational instability between these periods. 

Insulation may fail, while the carbon-steel pipe is absolutely in good shape, this calls an attention for continual maintenance of 
the insulator in other to protect the pipe from corroding and the best way to optimize instrumentation and process control is to 
reduce the cost of maintenance strategy. Thus, a single model which can predict insulation failure on a natural gas transmission 
pipeline at the exact spot along the pipeline where this failure will occur on a real time basis has been modeled. These showcase the 
accurate distance measures and ensure timely intervention on the affected area, initiates a maintenance shutdown and record “time 
between failure and time between repairs (TBF/TBR)” thereby preventing pipeline failures and reduction of operation operational 
instabilities. This model can be used to validate the readings of an Advanced Fiber-Optic Distributed Temperature Sensors (DTS) as 
well as Distributed Acoustic Sensors (DAS) in a gas transmission line. 

The complex study of instrumentation and process control, heat transfer, fluid mechanics and aerodynamics formed the basis at 
which this research was carried out. The LNG operational data which was subjected to Kalman filter and IQ-Rm PRO APIS software 
was also used to filter and generate data estimate and sensor reading analysis respectively. 

The natural gas pipeline operational data was filtered by kalman filter software and most of these data were used to model the 
position of the pipe, the flow velocity along the pipeline and exact spot where insulation failure will occur through sensors, 
Extendsim9 3D camera window, computerized display block which calculates the exact spot where insulation failed and Decision 
block known as indicator which uses logic as one of its tool on the gas arrival temperature to control its decision. 


Comprehensive Summary of Simulation Model development 
The following developed models were used as the programing code for the software: 


é 10 1 0 é 10 10Q, 06Q 
_— k))+——(rphV, )+-—(phV, )+—(phV. )+} —— —— at 
5 (ole+k))+—<(rahv, + (ohyV,)+<(ohv.)+| << (r9, + 2 +E 
av, (2% “| o) é (% ,) (& &) 
Cp + Tog a FT, . —| To) 7 ot +T,. —+ 
Or rood fr Oz or er r 00 Or 


10V, OV, Le Cty. Fm . CT, 1 OT gg 
——i4—* |/-V | —- +——¢ __-# 4 2 |_y | — 
ra(? 00 ve) (2%, (0% rood r © \r 00 
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1 OT 1 1). 1@¢,... Of 
ans) 2 sail pen -0 6: Zz 
ae of (r Tg)}+ Bc (2%olr +4} +2 3p + ae Energysource 


The model in equation (36) was derived from continuity, momentum and energy equations using the stated assumptions; and 


was written in C++ programing code in the simulation software which formed the basis for the 3-Dimensional internal flow: 


q= T, -T; 
7 a ua 
r. 
: +L, fh +L,“ +L, Bey 2 
QnhrL  2mk,L* 2nk,L* 2nk,L 2ahyrL 
T. -T, 
q= “ 


1 LL rl LL ri 3, 1 


QahrL * 2nk\L * 2nk,L 2ahr,L 


The model in equations (42a), (42b) was derived from conservation of energy equation and heat transfer models after the stated 
assumptions were made and was also written in C++ programing language code used to simulate the total heat loss for the Normal 
insulated/functioning system and that of the failed or faulty part of the system along the pipeline 


P=P,—puy 
_ ar 
Pe Ty=7 4 —m) 


These models in equation (45) and (48) were derived from conservation of momentum of the fluid element equations and 
conservation of energy equation after the stated assumptions were made. These also is written in C++ programming language Code 
used to simulate and predict the temperature/pressure profile for the Normal insulated/functioning system and that of the faulty or 
failed part of the system along the pipeline respectively. 


C,dT _ kAOT 


ot oror 


While the model in equation (51) was derived from the combination of conservation of momentum, energy equation and heat 
transfer after the stated assumptions has been made. The model is written in programing code to achieve model sensitivity on the 
flow problem. 

Due to the complicated nature of the models, Extendsim9 software is considered as the most suitable software to handle these 
problems numerically to achieve best possible results. 

However, the above equations is solved simultaneously with the models describing the device behaviour (compressor, maps, 
valve equations, set points, choked or free flow, equations for simulated failures, resistance equations, etc.),which are the most 
accurate. 


Methods used to validate this model 

The simulation software is made up of discrete rate, discrete event and continuous model blocks which have series of developer's 
tool to enable any researcher write or develop its model for any type of real life problem ranging from engineering, architecture, 
medicals, security, aviation, manufacturing, productions etc. This model is a combination of these three. 
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The tank block which represent source in discrete model is regarded as the refinery where raw natural gas flow through the 


gathering pipeline has an infinite capacity and was set to infinity since gas flows from the refinery are continuous. This means that 
natural gas beneath the earth crust has an infinite flow capacity and therefore its exploration is always continuous 

The flow of this natural gas is compressed using the compressor block tool and is being controlled by the valve block tool at the 
rate of 10000m?/day. The gas flowed through the pipeline where it is diverged by diverge block tool which is set to proportional at 
the diverge proportionality of 1:1 to diverge the flow at the rate of 5000m?/day to each segment and merge rate set to priority of 
outflows to merge the gas flows to different gathering pipelines down to the processing unit. At the processing plant, the valves are 
set to 3000m?/day to control the rest of the flows to the transmission pipelines. 

These merge and diverge block tools depicts its name in any model by diverging and merging of flows. Since tanks are used as 
source, it is also used as storage facility in Extendsim9 since it is known as residence block. The use of throw and catch block here is 
to minimize the number of blocks on this particular flow model. The good thing about this software is that it uses anything one can 
think of as well as defined within the acceptable range as block including ‘text’. 

The gas flow on the transmission pipeline has its flow rate at 3120m?/day as the flow along the transmission pipelines are set to 
sensing flows and priority of outflows respectively. The sensor is set at every 2km along the transmission pipeline where constraint 
variables or unstable variables of aerodynamics are read and recorded by the look up table block otherwise known as constraint 
variable tool, through this recordings, it sends the recorded result to the decision block otherwise known as the indicator to 
compare and analyze the readings of the temperatures, density etc. Extendsim9 understands the language of Computational fluid 
dynamics (CFD) and therefore uses true or false strings and attributes to differentiate the range at which heat transfers either by 
convection or conduction, and their temperature gradients. 

After the analysis, the decision block otherwise known as the model indicator is used to send feedback to both sensors and 
valves to initiate control of gas flow by closing the failed line arm and opening another arm of flow to diverge and merge flow to 
other pipes while automatically stopping gas flow to the affected area and also sends signal for a maintenance shutdown review. 
The computerized display function block is used to review the readings, calculate the failure distance from its current location and as 
well send feedback to the plotter and the indicator to declare the affected section/area ‘failed’ and therefore trigger shutdown repair 
while controlling the rest of the flow process to its final destination. With this, it sends feedback to the computerized display block 
to calculate and display the exact spot along the pipeline where the insulation failed. Every other block like interchange block, 
Resource item block, get, activity block, queue block are used for accuracy and accurate calculations to assist the decision block take 
its decisions. They also contribute to the successful flow of gas in other pipes along the pipeline. 

The results were displayed in Tables and Figures and properly discussed. 
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Figure 6 Normal/Field Operational Pressure versus Pipeline Distance 


Figure 6 illustrates the relationship between Normal/Operational pressure and pipeline distance upon the influence of dynamic 


characteristics in the flow line. Decrease in pressure in pressure was observed with increase in pipeline distance, owing it to frictional 
resistances due to body and surface forces in terms of pressure forces and other environmental factors from our findings. 
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Figure 7 Normal/Field Operational Temperature versus Pipeline Distance 


Figure 7 shows the NNPC operational data model in terms of temperature. From the graph, there is a decrease in temperature 
along the transmission line owing it to low viscosity in terms of temperature; thus, the initial and the arrival temperature and 
pressure is 51.172°C- 44.45°C. 
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Figure 8 Model Pressure versus Pipeline Distance (without failure) 
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Figure 9 Model Temperature versus Pipeline Distance (without failure) 


pase 6 1 


ARTICLE 


Figure 8 shows the relationship between modeled pressure and pipeline distance upon the influence of dynamic characteristics 


as the predictive model depicts the operational pressure. from our findings, it shows that the model was able to assert the readings 
of operational data as well as deduction of dynamic pressures at each kilometer along the transmission line. 

Figure 9 shows the relationship between the model temperature and pipeline distance upon the influence of dynamic velocity in 
the flow line as the predictive model depicts the operational temperature. Therefore this model can be used to predict the future 
behaviours of any pipeline system. 


Sensitivity Model and Normal Sensor readings Comparison 
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Figure 10 Model/Normal Operational Temperature verses pipeline distance 


Figure 10 shows the sensitivity comparisons for both model and operational sensor readings in terms of temperature. From the 
graph, we observed that the sensitivity model and sensitivity normal/operational temperature is almost the same. This shows the 
rate at which the model is sensitive to real-life Computational fluid dynamics (CFD) readings with calculated sensitivity rate (Sp) of 
0.922, which is 92.2%. This shows that the model can predict insulation failure along the transmission line with 0.982 accuracy i.e. 
92.2 degree accuracy. 
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Figure 11 Sensitivity graph of normal/operational and failure sensor readings versus time 
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Figure 11 shows the component sensitivity readings for both failure and normal sensor readings in terms of temperature and 


viscosity. From the graph, we observed that at the failed part of the system, the temperature readings increased with an increase in 
viscosity while the sensor readings of the normal system (without failure) maintains a normal flow progression owing it to body and 
surface increase. From our findings, it shows that the failed or faulty part of the flow system maintained a rapid increase towards the 
exit of the flow owing it to heat transfer by conduction from the surrounding to the pipe wall which transfers to fluid flow. Also the 
graph showed the variance between the operational/normal sensor readings and the failure sensor readings which means that the 
model is sensitive to both large and small scale variation on the pipeline reading which enable it to generate error and detect fault 
along the transmission line. 


a Failure Graph 
£ 70 
) 

F 60 
© 

2 50 
<= 40 

1S) 

@° 30 
& 20 
oO 

E 10 
fo} 0 
2 


1 3 5 7 9 11 13 15 17 19 21 23 25 27 29 31 33 35 37 39 
Pipeline Length Km 


Figure 12 Graph of Normal/Failed arrival temperature versus Pipeline distance 


The predictive model was used to predict the insulation failure along the transmission line using the arrival temperature as 
shown in Figure 12. This shows the arrival temperature for both normal (without failure) and the failed systems. From the graph; we 
observed that the normal (without failure) system maintained its arrival temperature level towards the exit of the flow while the 
failed or faulty part of the system increased its arrival temperature exponentially from 12.5km and maintained its increase towards 
the exit flow. From our findings, the increase found at the failed or faulty part of the system is as a result of erosion or pitting 
(leakage) found along the transmission line which made the fluid temperature depicts the surrounding temperature hence heat is 
transferred and because the rest of the pipe wall is well insulated, it maintained an increased temperature towards the exit of the 
flow. 
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Figure 13 Graph Illustrating the Exact Location of the Failure 
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Once more, if the temperature of the surrounding is below the normal (without failure) system temperature, the failed part will 


also adapt to the surrounding temperature. For instance; in arctic region where we have a more freezing temperature, since gases 
are cryogenic in nature, tends to freeze if it adopt the surrounding temperature and this in particular leads to fogging or hydrate 
formation thereby causing pipe blockage along the transmission line. However, if the gas temperature adapts to a high surrounding 
temperature, with an absolute reoccurrences, it may lead to explosion, this calls for frequent checks on our gas transmission line in 
order to maintain the pipeline as at when necessary to avoid ‘guess work’ during inspection. 

Figure 13 demonstrates failed area and pipe location against time shows the exact area of the failed insulator/leakage along the 
pipeline. From the computerized displayed function block (tool), it was shown that the calculated distance and the exact spot of 
failure along the pipeline using monitor indicating sensor which reads 12500m (12.5Km) due east from the simulated metering 
station (office). Here the exact spot where insulation failed along the pipeline is 12500m due east. 


Plate 2 Extendsim9 3D Camera Model (Subsea Visual View for Exact Spot of Insulation Failure) 


Plate 2 shows the 3D visual view of the failed insulator along the natural gas transmission for offshore pipelines using 
Extendsim9 3D phase from our observations, the dark part of the plate shows the well-insulated area while the spot-light areas are 
the failed spots or eroded part of the pipeline system. The area is covered with erosion of insulator and pitting/leaks. 
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Plate 3 Extendsim9 3D Camera Model (Underground/Onshore Visual View for Exact Spot of Insulation Failure) 
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Plate 3 shows the 3D visual view of the failed insulator along the natural gas transmission for onshore pipelines using 
Extendsim9 3D phase from our observations, the dark part of the plate shows the well-insulated area while the spot-light areas are 
the failed spots or eroded part of the pipeline system. The area is covered with erosion of insulator and pitting/leaks. 
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Sensor 1 reading 


Plate 4 3-D shows the distance between the eroded spot and pitting along the failed distance 


Plate 4 shows the cross section of the failed area and the distance between pits using the software 3D camera view from sensor 
readings. From our observations, the model was able to detect the spots where insulation failed along the pipeline: 


Sensor % Feadily 


Plate 5 Sensor one and sensor two position reading 


Plate 5 illustrates the model 3-D camera view; we observed that gas produces more noise as it approaches the failed area more 


especially within the erosion area as it expand from its control volume and decrease the noise level at the failed spot (leak or pit), 
these act as a warning signal during operation. 


4. CONCLUSION 


Generally, in gas flow measurement, the density of the gas changes as the pressure and temperature change. This change in density 
can affect the accuracy of the measured flow rate if it is uncompensated. When gas flow through pipes, it tends to experience 
pressure drop which occur as a result of friction resistance along the natural gas transmission line. Since gases are compressed, its 
density is a function of absolute temperature and absolute pressure, the change in density is proportional to the change in pressure 
or temperature. However, the study of instrumentation and process control, heat transfer and fluid mechanics formed the basis of 
this model. In this dissertation, a model that can predict insulation failures along natural gas transmission pipelines was developed, 
validated using sensitivity analysis which shows that the model can predict insulation failure on a real time basis since the rate of the 


Pane 65 


ARTICLE 


model sensitivity is 0.982. From the analysis, it showed that the arrival temperature increased exponentially across the failed area 


due to the erosion and or leakages as the gas adopts the surrounding temperature and maintained a temperature increase towards 
the exit flow. This is due to the fact that the rest of the pipe walls are well insulated. In other words the model can detect failure 
using the arrival temperature. The display tools also displayed 12.500Km due east from its location as the exact spot of the failed 
insulator and a maintenance shutdown from the model was used to record time between “failures and time between repair”. Also 
from the analysis, it showed that this model can read and interpret distributed temperature sensing devices or distributed acoustic 
sensing device along the pipeline. The result showed that this model can predict insulation failure along natural gas transmission 
pipeline on a real-time basis by analyzing the flow rate and temperature of the flowing gas along the line. This will solve the late 
detection of insulation failure problems for both insulation failures with/without pipe deterioration. 


Funding: This research received no external funding. 


Conflicts of Interest: The authors declare no conflict of interest. 


REFERENCE 

1. Chetounani, Y. (2008). Design of a Multi-model Observer- 
Based Estimator for Fault Detection and Isolation (FDI) 
Strategy Application to a Chemical Plant. 10-48 

2. Dasjalais, O. A. and Zarr, R.R. (2002) Insulation materials 
testing and applications, New York ASTM_ international. 
Wikipedia, free encyclopedia 

3. Gajanan, D., Preeti, B., Rupesh, D. and Saurabh, P. 
(2017).Thermal Insulation Material Tools for Engineering 
Conservation: |ICAR-National Diary Research Institute SRS, 
Bangalore; Indi, 1-40. 

4. Guo and Ghalambor, (2002).Natural Gas Engineering 
Handbook. New York: Gulf Publishing Company, 1-30 

5. Lawal, M. O. (1989).Pipeline Development in the Nigerian 
Petroleum Institute. Geography, 71(1): 60-62. 

6. Li, C., Jia, W., and Wu, X. (2011) Modeling and Simulation for 
Steady-State and Transient Pipe Flow of Condensate Gas: 
School of Petroleum Engineering, Southwest Petroleum 
University, China: 66-80 

7. Ogwu, F. A. (2011). Challenges of Oil and Gas Pipeline 
Network and the Role of Physical Planners in Nigeria, 
FORUME journal: (10): 41-51 

8. Timur, C. (2007) Pipeline Leak Detection Techniques: Andrei 
Saguna University, Constana Romania, Computer Science. 
Constana Romania: Anale Seria Informatica. 

9. Yung-Sang, K. (2018). Engineering Failure Analysis: School of 
Advanced Materials Science and Engineering: 
Sungkyunkwan University Republic of Korea. 

10. Yung-Sang, K. and Jung-Gu, K. (2017).Failure Analysis of a 
Thermally Insulated Pipeline in a District Heating System. 4- 
44 

11. Zarr, R. R. (1989). Steel Pipe: A Guide for Design and 
Installation. New York: American Water Works Association. 


pane | 66 


